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Abstract 

The dynamics of a cascaded system that consists of two atom-cavity subsystems is studied by using the quantum trajectory 
method. Unwanted losses are included, such as photon absorption and scattering by the cavity mirrors and spontaneous 
emission of the atoms. Considering an initially excited two-level atom in the source subsystem, analytical solutions are 
obtained. The entanglement evolution is studied for the two atoms and for the two intracavity fields. 
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1. Introduction 

The concept of entanglement has been of great in- 
terest since the early days of quantum mechanics [1], 
and it has become of central importance in a vari- 
ety of discussions on the fundamental aspects of the 
theory [2,3]. Unexpectedly, entanglement has also 
been demonstrated in classical Brownian motion [4]. 
Nowadays entanglement is receiving new attention 
in the rapidly developing fields of quantum informa- 
tion, quantum computation and quantum technology; 
for reviews, see [5,6]. In the context of entanglement 
preparation between atoms at separate nodes, a vari- 
ety of schemes have been proposed, for example, by 
measuring the superpositions of light fields released 
from separate atomic samples or by measuring a probe 
light field that has interacted in a prescribed way with 
different samples. Due to the indistinguishability in 
the measurement and conditioned on the results of the 
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measurements, the atomic system is projected onto 
an entangled state [7,8,9]. An unconditional prepara- 
tion of entanglement has also been analyzed in the 
case of a cascaded system. This has been discussed 
for two distantly separated atoms [10,11], as well as 
for separate atomic ensembles [12]. More recently, the 
entanglement evolution for a Raman-driven cascaded 
system has also been analyzed [13] . 

In the spirit of these previous achievements, in the 
present contribution we will consider a cascaded open 
quantum system. We study the dynamics of a system 
that consists of two atom-cavity sub-systems A and B. 
The quantum source A emits a photon and the second 
quantum subsystem B reacts on the emitted photon. 
Unwanted losses are included, such as photon absorp- 
tion and scattering by the cavity mirrors and sponta- 
neous emission of the atoms. Considering an initially 
excited two-level atom in the source subsystem, ana- 
lytical solutions are obtained. Subsequently, the entan- 
glement evolution between the two atoms, as well as 
between the two intracavity fields, is studied by using 
the concurrence. 
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The paper is organized as follows. In Sec. 2 the mas- 
ter equation describing the dynamics of the cascaded 
system is introduced, and the problem is solved analyti- 
cally by using the quantum trajectory method. In Sec. 3 
the entanglement evolution between the two atoms, or 
between the two intracavity fields, is analyzed. Finally, 
some concluding remarks are given in Sec. 4. 



2. Cascaded system dynamics 

In this section we analyze the dynamics of the sys- 
tem under study. The cascaded open quantum system 
consists of two atom-cavity subsystems A and B, where 
the source subsystem A is cascaded with the target 
subsystem B, cf. Fig. 1. The cavities have three per- 
fectly reflecting mirrors and one mirror with transmis- 
sion coefficient T -C 1. In the two subsystems A and B 
we consider a two-lovcl atomic transition of frequency 
LOk (related to the atomic energy eigenstates \lk) and 
|Ofc)) coupled to a cavity mode of frequency w^, where 
k = a,b denotes the subsystem. The cavity mode is do- 
tuned by Afe from the two-level atomic transition fre- 
quency, ujk = Wfc -h Afc , and is damped by losses through 
the partially transmitting cavity mirrors. In addition 
to the wanted outcoupling of the field, the photon can 
be spontaneously emitted out the side of the cavity into 
modes other than the one which is preferentially cou- 
pled to the resonator. Moreover, the photon may be 
absorbed or scattered by the cavity mirrors. It has been 
shown that unwanted losses can have significant effects 
on the dynamical evolution of a quantum system and 
cannot be a priori neglected, see, e.g., Refs. [14,15]. 
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Fig. 1. The cascaded open system consisting of two atom- 
cavity subsystems A and B. 



To describe the dynamics of the system we will use a 

master equation formalism. This leads to the following 
master equation for the reduced density operator p{t) 
of the system: 



dm 

dt 



JimJi-^JlJim 



(1) 



The Hamiltonian is given by 



H = Ha + Hb + ih^^^ [e-*na^ - e'^b^a) , (2) 

where Ha and Hb describe the atom-cavity interaction 
in the two subsystems A and B, respectively, and, in 
the rotating-wave approximation, are given by 

Ha = hga [aAio + a^Aoi) -l- /jA^iii , (3) 

and 



Hb = hgb [bBw + b^Boij + hAbBn . (4) 

The third term in Eq. (2) describes the coupling be- 
tween the two cavities [16,17]. In these expressions, a 
(a^) is the annihilation (creation) operator for the cav- 
ity field A and similarly b (6^) for the cavity field B. 
We have also defined Aij = \ia){ja\ = 0, 1), and 
Bij = |j6){ib| = 0, 1). In addition, is the atom- 
cavity coupling constant and nk the cavity bandwidth. 
The phase 4> is related to the phase change upon re- 
flection from the source output mirror, and/or to the 
retardation of the source due to the spatial separation 
between the source and the target, cf. [18]. 

The jump operators ■/» entering in Eq. (1) are defined 

Ji = sfi^a + y/Kbe~'''^b , (5) 
which describes photon emission by the cavities; 



J2 



J3 



(6) 



are associated with photon absorption or scattering by 
the cavity mirrors; and 



Ja 



(7) 



are related to a photon spontaneously emitted by the 
atoms. Here n'k and Vk are the cavity mirrors' absorp- 
tion (or scattering) rate and the spontaneous emission 
rate of the two-level atom, respectively. Note that the 
operator Ji contains the superposition of the two fields 
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radiated by the two cavities due to the fact that the ra- 
diated photon cannot be associated with photon emis- 
sion from either A oi B separately. 

In the following we will identify, for notational con- 
venience, the state \a) with the state |1, 0, 0, 0), which 
denotes the atom A in the state |la), the cavity A in 
the vacuum state, the atom B in the state \0b) , and 
the cavity B in the vacuum state. In the state \b} = 
|0, 1, 0, 0) the atom A is in the state |0a), and the cavity 
A is in the one-photon Fock state. Similarly, we define 
\c) = 10,0,1,0), \d) = |0, 0,0,1), and |e) = |0. 0,0,0). 
The state \a) will be considered as the initial state 
of the system. It follows that the Hilbert space that 
describes the cascaded system under study is, in our 
model, sparmed by the five state vectors |a), |6), |c), 
\d), and |e). 

To evaluate the time evolution of the system we use 
a quantum trajectory approach [18,19,20]. Let us con- 
sider the system prepared at time to = in the state 
I a). To determine the state vector of the system at a 
later time t, assuming that no jump has occurred be- 
tween time to and t, we have to solve the nonunitary 
Schrodinger equation 

*S^IV5no(t))=J?'|V;no(t)), (8) 

where H' is the non-Hermitian Hamiltonian given by 



H' =H-'^J2jlJi=HA + HB- ih^^a^a 

i=l 

+ ^b^b + ^iii + y Sii + V^e^'^&ta) , (9) 



where we have defined 

Ka = K,a + l^a : Kb — l^b + l^'b ■ 



(10) 



If no jump has occurred between time to and t, the 
system evolves via Eq. (8) into the unnormalized state 

|V;„o(t)) = a{t)\a) + /3(t)|6) + ^{t)\c) + 5(t)|d) . (11) 

In this case the conditioned density operator for the 
atom-cavity system is given by 



Pno(t) 



|Vno(t))('(/'no(t) 



(12) 



(Vno(t)IVno(t)) ' 

where we indicate with "no" the fact that no jump has 
occurred between time to and t. 

The evolution governed by the nonunitary Schrodinger 
equation (8) is randomly interrupted by one of the 



five kinds of jumps J,, cf. Eqs. (5)-(7). If a jump has 

occurred at time tj, tj € (to,t], the wave vector is 
found collapsed into the state | e) due to the action of 
one of the jump operators 



J^ IV'no(tj)) -> |e) (i : 



,5). 



(13) 



In the problem under study we may have only one 
jump. Once the system collapses into the state |e), the 
nonunitary Schrodinger equation (8) lets it remain un- 
changed. In this case the conditioned density operator 
at time t is given by 



Pyes(t) = |e)(e| , 



(14) 



where we indicate with "yes" the fact that a jump has 
occurred. 

In the quantum trajectory method, the density op- 
erator p(t) is obtained by performing an ensemble av- 
erage over the different conditioned density operators 
at time t, yielding the statistical mixture 



p(t) = Pno(t)pno(t) +Pycs(i)/5yes(t) . 



(15) 



Herepno(t) andpyes(t) are the probability that between 
the initial time to and time t no jump and one jump 
has occurred, respectively, where Pno(t) -l-Pyes(t) = 1. 

To evaluate Pno(t) we use the method of the delay 
function [20]. This yields the probability Pno(t) as the 
square of the norm of the unnormalized state vector: 



Pno(t) = II |'0no(t)) ||^= ('0no(t)|Vno(t)) 

= |a(t)|^ + |/3(t)|^ + |7(t)|=^ + |5(t)|^ 



(16) 



From Eqs. (15) and (16) one obtains for the density 
operator p{t) the expression 



p(t) = |V;no(t))(V'no(t)| + |e(t)|'|e)(e| 



(17) 



where we have defined |e(t)p = 1— Pno(t). The physical 
meaning of |a(t)|2, \f5{t)\^ , \^{t)\^, \5{t)\\ and |e(t)|2 is 
clear. They represent the probability that at time t the 
system can be found either in \a), \b), \c), \d), or |e). 

In order to determine Q:(t), /3(t), 7(t), and (5(t), we 
have to solve the nonunitary Schrodinger equation, cf. 
Eqs. (8) and (9). This leads us to consider the following 
inhomogeneous system of differential equations, similar 
to the one in Ref. [13], 
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a{t) = -i (Aa - iTa/2) ait) - igaP{t) , 

m = -igaa{t) - {Ka/2)m , 

jit) = -i (Ab - iTb/2) 7(t) - igbSit) , 

Sit) = -igbiit) - iKi/2)Sit) - ./J^e^^m . 

(18) 

The differential equations for ait) and /9(t) can be 
solved independently from those for jit) and Sit). For 
the initial conditions a(0) = l and /3(0) = 0, and defin- 
ing 



fife = ^ ^-4gl-iKk -i^) - (^Afe - 



(19) 



we can write the solutions as 
"if„/2-i(Aa-ir„/2) 



ait) = 



sinh 



^at 



+ cosh 



0.at 



-[(ifa+r„)/4+iA„/2]t 



Pit) = - ^ sinh f g-[(K,+r,)/4+iA„/2]* 



(20) 



Inserting now in the inhomogeneous pair of differential 
equations for fit) and Sit) the solution obtained for 
/3(t), wo can determine the solutions for 7(i) and Sit). 
For the initial conditions 7(0) = and 5(0) = 0, we can 
write the solutions as 



7(i) = gb {f+it) [g- it)+h+it)]-f- it) [g+ it) + h- it)]} , 
Sit) = i 



Ki-Ti .At fii, 

— *-+T 

Kj-Tb .At fib 
4 * 2 2 



f-m+{t)+h-it)] 

f+it)[g-it)+h+it)],i21) 



where we have defined, for notational convenience, 



^ ga^yUgKbe' ^[-(K|,+rb)/4-iAfc/2±fifc/2]« ^ ^22) 

[(n„±n6)/2-T-iA]t _ -, 
9± (t) = ^/o ^ . (23) 



(Oa ± fi6)/2 - T - iA 



and 



g-[(r!a±fi!,)/2+T+<A]t _ 

= (n„±nt)/2 + T + iA ' 

where T = (jr,-7^6+r„-ri,)/4 and A=(Aa-A6)/2. 
In the case of equal parameters for the two subsystems 
A and B, the solutions (21) simplify as 
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Fig. 2. The probabilities |Q(t)P (dotted line), |/3(t)P (dot- 
dashed line), |7(t)P (full line), and |<5(t)p (dashed line) are 
shown for the case of equal parameters for the two subsys- 
tems A and B, where g/K = 5, k/K = 0.9, A/K = 0.1, 
T/K = 0.2. 



-rit)- 



Kg^e 



j-(if+r)/4-iA/2+n/2]t 



K-r A n 

^-*T+2 



V J-(K+T)/4-iA/2-n/2]t ll^ge 



Ot- 1 



K-V .A 
^-*2 



[e-"*+m-l] e[-(^+r)/4-^A/2+n/2]t^ (35) 



where wo have used limj:^o{[oxp(±a;f) — l]/a;} = ±f, 
and defined k = Ka = Kb, K = Ka = Kb, A = Aa = At, 
T = Va = Tb, g = ga=gb, and Q, = Q.a=flb. In Fig. 2 we 
show the functions |a(^)|^ l/3(t)P, |7(t)P, and \5it)\^ , 
i.e. the occupation probabilities of the state | a) , 1 6) , | c) , 
and |d), respectively, for the case of equal parameters 
for the two subsystems A and B, mthg/ K = h, n/ K = 
0.9, A/K = 0.1, and T/K = 0.2. For these functions 
the phase factor e"* does not play any role. 



3. Entanglement evolution 

In the system under study, the two atoms constitute 
a pair of qubits. An appropriate measure of the entan- 
glement for a two qubits system, often considered in 
the context of quantum information theory, is the con- 
currence [21]. Given the density matrix p for such a 
system, the concurrence is defined as 
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Fig. 3. The concurrence C[pat(t)] between the two atoms is 
shown for the case of equal parameters for the two subsys- 
tems A and B, where g/K = 5 and A/A' = 0.1. Moreover, 
k/K = 0.9 and T/K = 0.2 (sohd hne), k/K = 0.9 and 
r/K = (dotted line), k/K = 1 and T/K = (dashed 
line) . 



C{p) 



{o, 



0^2 



/X. 



I}, (26) 



where Ai > A2 > A3 > A4 arc the eigenvalues of the 
matrix p = p{<yy®<Jy)p*{(jy ®<jy). Here ay is the Pauli 
spin matrix and complex conjugation is denoted by 
an asterisk. The concurrence varies in the range [0, 1], 
where the values and 1 represent separable states and 
maximally entangled states, respectively. 

To derive an expression for the concurrence between 
the two atoms, let us consider the density operator that 
describes the system. It is obtained from the density 
operator p(t), Eq. (17), by tracing over the intracavity 
field states for the two subsystems, pat [t) = Trcav [p{t)] , 
and is given by 

Pat(t) = |a(t)ni,0)(l,0| + |7(t)nO, 1)(0,1| 

+ a(i)7* (i) 1 1, 0) (0, 1| + a* (t)7(t) |0, 1) (1, 0| 

+ {l-[Ki)|^ + |7(t)n}|0,0)(0,0|. (27) 

Considering the 4x4 density matrix pat it) , it is easy to 
show that the concurrence C[pax(t)] is, using Eq. (26), 
given by 

C[pat(t)] = 2|a(t)||7(t)| . (28) 
To analyze the time dependence of the concurrence 
between the two atoms, let us consider the case of equal 
parameters for the two subsystems A and B. Inserting 
the analytical solutions (20) and (25) into Eq. (28), we 
show in Fig. 3 the function C[pat (t)] for the parameters 
g/K = 5, A/A' = 0.1, and for different values of n/K 
and V /K. Since the concurrence contains only absolute 
values, the phase factor e"'' does not play any role here. 



Fig. 4. Comparison between the concurrences for the two 
atoms (solid line) and for the two intracavity fields (dashed 
line). The parameters for the two subsystems are equal, 
with g/K = 5, A/ZsT = 0.1, k/K = 0.9 and T/K = 0.2. 

From Fig. 3 one can clearly see that the initially 
disentangled atoms become entangled. In particular, a 
maximum value for C[pat(t)] is found for i ~ 1.88/Jsr, 
where, for the shown cases, C[pat(f)] — 0.74, 0.66, and 
0.55. iFxoxa the entanglement between the two atoms, 
as shown in Fig. 3, it is clearly seen that it sensitively 
depends on unwanted losses. Already the effects due 
to the absorption or scattering by the cavity mirrors 
alone are not negligible. Moreover, the inclusion of the 
spontaneous emission further reduces the degree of the 
entanglement between the two atoms. For example, the 
relative variation of the concurrence is approximately 
10% between the ideal case it/K — 1 and T/K = (no 
absorption or scattering and no spontaneous emission) 
and the case k/K = 0.9 and V/K = 0, considering 
the peak at t. The concurrence decreases even further 
if spontaneous emissions are included. In this case for 
k/K — 0.9 and T/K = 0.2 the peak of the concurrence 
at i is reduced by approximately 25% compared with 
the ideal lossless case. These values clearly tell us that 
the effects of unwanted losses cannot be neglected in 
general, when considering the entanglement in realis- 
tic quantum systems. Of course, for Kt ^ 1, the two 
atoms become again discntaugled due to the omission, 
sooner or later, of a photon into one of the five decay 
channels. The release of a photon into the environment 
destroys any entanglement, projecting the two-atom 
subsystem onto the separable state |0, 0). 

Finally, we note that the concurrence between the 
two intracavity fields can be obtained as well. Let us 
consider the density operator that describes the sys- 
tem of the two intracavity fields A and B, obtained 
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from the density operator p(t), cf. Eq. (17), by tracing 

over the atomic states of the two subsystems, pcav (t) = 
Trat [p{t)]- Considering now the 4x4 density matrix 
pcav(t) in the two intracavity-field Fock basis, the con- 
currence C[pcav(t)] is given by 
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c[po.At)] = '2\m\m\ ■ (29) 

A comparison between the concurrence for the two 
atoms and the concurrence for the two intracavity fields 
is shown in Fig. 4. When the concurrence between the 
two atoms reaches a majdmum, the concurrence be- 
tween the two intracavity fields is approximately zero, 
and vice versa. This is related to the fact that the ex- 
citation energy is transferred between the atoms and 
the intracavity fields. When the atoms are unexcited, 
their state is separable and the intracavity fields are 
entangled. As a function of time, the entanglement is 
thus exchanged between the two atoms and the two 
intracavity fields. 



4. Summary 

The dynamics of a cascaded system that consists of 
two atom-cavity subsystems has been analyzed. Un- 
wanted losses have been included, such as photon ab- 
sorption and scattering by the cavity mirrors and spon- 
taneous emission of the atoms. The evolution of the 
open quantum system under study has been described 
by means of a master equation. Considering an initially 
excited two- level atom in the source subsystem, ana- 
lytical solutions for the dynamics of the system have 
been obtained. The entanglement evolution between 
the two atoms, constituting a two-qubit system, has 
been studied by using the concurrence. A similar anal- 
ysis has been performed for the two intracavity fields. 

The dynamical evolution of the system shows that 
the two initially disentangled qubits reach states of 
significant entanglement. It has been shown that the 
entanglement generated between the two atoms sen- 
sitively diminishes due to the presence of unwanted 
losses, which cannot be neglected in realistic quantum 
systems. We have also shown that the entanglement is 
periodically exchanged between the two atoms and the 
two intracavity fields. 
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